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Abstract 

We discuss some new invariants of quark mixing and show their usefulness with a simple 
example. We also present some other new tools for analyzing quark mixing. 
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1 Invariants of quark mixing and their applications 



In recent years the increasingly accurate data on quark mixing has stimulated interest 
in possible predictions concerning quark mass matrices. The parameters of the Cabibbo- 
Kobayashi-Maskawa mixing matrix V are now known reasonably well This determina- 
tion has been made possible partly by the finding that there are only three generations of 
usual standard- model fermions (with corresponding light or massless neutrinos). Since the 
diagonalization of the quark matrices in the up and down sectors determines V, one can 
work back from the knowledge of V to put constraints on the possible forms of (original, 
nondiagonal) quark mass matrices. However, the data on quark mixing determines these 
mass matrices only up to an arbitrary unitary similarity transformation. This is a result of 
the fact that if the up and down quark mass matrices, M„ and M^, are both acted on by the 
same unitary operator Uq according to 

Mu,d ^ UoMu,dU^o (1) 

then the mixing matrix V remains unchanged. There have been many attempts to study 
specific assumed forms for quark matrices. While this is worthwhile, it is desirable to express 
the constraints from data on V on the quark mass matrices in an invariant form. In Refs. 
0, 1^ certain invariant functions of the quark mass matrices Ipg were introduced, which are 
expressed in terms of the quark masses squared and the \Vij\: 

I,, = Tr{Hl HI) = Y: {mt^f^imff^ \V.,\' (2) 

where Hq = M^M^, and mj"'' and m\^^ are the masses of quarks in the "up" and "down" 
charge sectors. 

In this paper we introduce some new invariants of the quark mass matrices (with respect 
to the transformation (|T|)) which can be expressed in terms of the measurable quantities only. 
These new invariants help one simplify the algebraic expressions which relate the elements 
of the quark mass matrices to the data. This important advantage allows for some new uses 
of the invariants which we will illustrate with some examples. 
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We introduce the following new invariants of the transformation (|l]): 



K„{a,f3) = det{aHl + (3HI) (3) 

where p,q,a, P 7^ 0. 

The hermitian matrices and can be diagonalized by a unitary similarity transfor- 
mation: 

where Dq = diag{{m^i^y , (m^^y) are the diagonal matrices of the quark masses 

squared. 

The mixing matrix V can be written then as: 

V = U^Ul (5) 
In order to find an expression for Kpq in terms of Uij we will need the following 
Theorem 

If A and B are two 3x3 matrices such that det{A) 7^ and det{B) ^ then the following 
relation holds: 

det{A + B) = det{A) + det{B) + det{A) Tr{A-^B) + det{B) Tr{AB-^) (6) 



Proof: 

We denote the elements of matrices A and B by Aij and Bij correspondingly. Their 
co-factors (which are equal to the corresponding minors, up to sign) will be written as Aij 
and Bij. Then each determinant may be decomposed in a sum (Laplace expansion): 

det{A) = ^ Aj-^ij = X! ^ij-^ij 

i j 

det{B) = '^BijBij = Bij Bij 

i j 
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By definition, the determinant of a 3 x 3 matrix is a sum of 3! = 6 terms: 

det{A + B) = J2{-lY{A,k, + B,kJ{A2k, + B2k,){A3k, + B^^,) (7) 

where r is the sign of the permutation (^^ \^ \^). 

The terms in the sum which contain only the elements of A can be arranged as det{A). 
Similarly, the terms containing only i?'s give det{B). The terms containing one element of 
A multiplied by two elements of B, or visa versa, can be rewritten as: 

+ (8) 

We can now use an identity: 
to rewrite (|) as: 

Y^iA.Bij + BijAij) = det{B) Y.Aij{B-^)ji + det{A) Y.{A-^),iBij = 

i,j ij ij 



det{A) Tr{A-^B) + det{B) Tr{AB-^) 

Altogether we get 

det{A + B) = det{A) + det{B) + det{A) Tr{A-^B) + det{B) Tr{AB-^) 

which is the statement of the theorem (|^). This completes the proof. 

Theorem @ may be easily generalized to the case of 2 x 2 matrices, in which case the 
last two terms in (^) are equal and correspond to a redundant counting of the same terms 
in a sum similar to (|^). Thus for the 2x2 matrices we get: 

det{A + B) = det{A) + det{B) + det{A) Tr{A-^B) = 
det{A) + det{B) + det{B) Tr{AB-^) 

The immediate consequence of equations @ and (H) is the following relation: 
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Kpg{a,(3) = det{aHP + f3H'^) = 
«3 (ml^^m^^^mi"^^ [ 1 + (/?/«) ^ [{mf^ fy{m't^)^P] U,, ] + (10) 

We also notice that 



K„{a, ±(3) = a3(mS"U5"^mi"))2^(l ± (/3/a)/(-rt,) 

±/?3(mS'')m('')mf )2''(1 ± («//?)/,(_,)) (11) 

Any four independent invariants from the set {Ipq, Kpg} contain all the physical informa- 
tion about the CKM matrix. 

If the mass matrices are assumed to be hermitian, one can introduce a similar set of 
invariants: 

i„ = Tr{MlMl) 

kpg{a, f3) = det{aMP + (3Mj) ^ ' 

The formulae, similar to (@), ([T0|) and (pT]), will also hold for J, K: 

4 = Tr(M„^ Ml) = Y: {mf'^rimfr \V.,\' (13) 

ij 



kpg{a,(3) = det{aMl + (3Ml) = 
a' (mS"^m^"U^"^)P [ 1 + {(3/a) ^ [(mf )7(mf ))^'] U,, ] + (14) 

f3' {mf^mf^mfy [ 1 + {a/ (3) ^ [{mfflimf^f] U,, ] 

The odd powers of the mass eigenvalues may appear in some of the / and ii'-type in- 
variants. As usual, the signs of the fermion masses are ambiguous. However, for a given 
model, different choices of signs will in general result in different predictions for the CKM 
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matrix. This is because in general the mass matrices do not commute with all of the diagonal 
matrices of the form diag{±l, ±1, ±1), where the + and — signs are chosen arbitrarily but 
so as to not get the plus or minus identity. We discussed the effect of the sign choices on 
the CKM mixing parameters on the example of some particular model 0. Once the choice 
of signs for the fermion eigenvalues (determined in practice by the best fit to the data) is 
made, there is no further sign-related ambiguity in the / and i^'-type invariants. 

Now we would like to illustrate the usefulness of the invariants discussed above. As 
an example, we will take the model proposed in ||^. This model gives predictions for the 
low energy data on fermion masses and mixing which are in reasonable agreement with 
experiment. This model is formulated in the context of an SO(IO) supersymmetric grand 
unified theory. The model has the following Yukawa matrices at the GUT scale: 






A-u 





A-u 


Bu 














Mu=\ Au Bu (15) 



( ° 
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Bd 
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-SBdC^' 


-SBd 





-Wd 


Cd 



(16) 



Me = Ade-'"^ -SBdC^' -SBd (17) 
V -3Bd Cd J 

The low energy effective theory below GUT thresholds is assumed to be the Minimal 
Supersymmetric Standard Model (MSSM). After the renormalization effects are taken into 
account, the model agrees with the data on fermion mixing for a broad range of the t-quark 
mass, rrit = 150... 190 GeV. In these fits the phase 6 is relatively unimportant and can be 
taken to be zero. 

The values of the parameters \Aq\ and Cq {q = u,d), in ([T5|) and (|l^) are simply related 



to the masses of quarks. On the contrary, the phase 0, which must have a nonzero value in 
order for the model to agree with experiment, depends on both the masses of quarks and 
their mixing parameters. It is a common procedure to diagonalize the matrices Mu and Md 
numerically (or by means of the Taylor series expansion in terms of the quark mass ratios) 



and compare the corresponding CKM matrix to the data. On the other hand, the method 
of invariants discussed above offers a simpler and more elegant solution: the phase can be 
expressed analytically in terms of the measurable quantities only. Because the mass matrices 
in ([T5|) and (|1^) are hermitian for ^ = 0, one can use a K-tjpe invariant. On one hand, the 
allowed range for the value of 1) is known in terms of the measurable quantities from 

(0) . On the other hand, 

ku{l,l) = det{Mu + Md) = 

= \A^ + A,\ (Cu + C,)= (18) 
= {2\A^Ad\cos{^) + + + Cd) 

And therefore 

iM I - _ (19) 

l^ul \^d\ 

The right-hand side of (|19D is known in terms of the quark masses and the CKM mixing 
parameters. Therefore, one can use the relation (^) to evaluate the phase without having 
to explicitly diagonalize the mass matrices. 

For this model and others, it is important to have a general methodology to determine 

(1) how many unremovable phases there are in quark mass matrices and (2) which elements 
of these matrices can be rephased to be real. We have presented a general method to answer 
both of these questions |^. The analogous questions for lepton mass matrices are more 
complicated, but we have also given a general answer to them 0. 

In summary, we have introduced and studied some new invariants of the quark mass 
matrices which are model- and weak basis-independent. The identities (|10D and (0), as well 
as (|r^) and ([T^), involving these invariants, provide important constraints on the possible 
forms of quark mass matrices and since they directly relate the elements of these 
matrices to the measurable parameters |Vijp and quark masses, and thereby enable one to 
avoid the explicit calculation of the eigenvectors of and Ha- 
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2 A method to generate families of viable quark mass 
matrices 



In this section we would like to present some useful mathematical tools, first proposed in [0, 
for the study of quark mass matrices and the connection with quark mixing. Our method 
was recently applied in p to generate the families of acceptable solutions for the fermion 
mass matrices. 

The idea of the method is to utilize the well-known experimental fact that the CKM 
matrix is close to the 3x3 identity matrix: \Vij\ ~ 6ij. We make use of this fact by writing 
V as 



V = e'""" (20) 

where H is some hermitian matrix and a is a real number. One may choose H to have its 
dominant (largest, in absolute value) eigenvalue to be 1. Then for a consistent with the data 
we find that |a| ~ 0.3. We can now expand V in the powers of a: 

V = l + iaH - \a^H'^ + ... + —AiaHf + ... (21) 

It was shown in Ref. that for any practical purposes it is sufficient to consider the 
first and the second order in a to match the precision to which the relevant quantities (CKM 
parameters, quark masses, etc.) are known. 

The matrix iJ, corresponding to a given matrix V with distinct eigenvalues fj, can be 
easily computed using the Sylvester's theorem: 

laH = }^ ln{vk)^ — J — 22 

Then the usual system of equations: 



Uu,LMuUlji = diag{mu,mc,mt) = Du 
< Ud,LMdUl ^ = diag{md, rus, rrib) = Dd (23) 
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is satisfied for tfie family of solutions: 

' M„ = UlDuUu = 
Du + iax[D^, H] - la''x^[[D^, H],H] + ... 

(24) 

Md = U^DdUd = 

+ ^a{x - l)[Da; H] - \a\x - lf[[Dd. H],H] + ... 

depending on some arbitrary parameter x. (It is assumed that as well as |1 — a;|, is 
sufficiently small to preserve the convergence of the series.) 

To summarize, in this presentation we have discussed two different approaches to ana- 
lyzing quark mixing. First, we introduced some new invariants of quark mixing and showed 
the usefulness of our method of invariants on a simple example. Then we have also discussed 
some new mathematical tools which allow one to generate the mass matrices consistent with 
the data. 

The author would like to thank Professor Robert Shrock for many helpful discussions 
and comments. 
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